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Abstract 

The quasi-local energy conservation law is derived from the vacuum Einstein's 
equations on the timelike boundary surface in the canonical (2,2)-formalism 
of general relativity. The quasi-local energy and energy flux integral agree 
with the standard results in the asymptotically flat limit and in spherically 
symmetric spacetimes. 
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In general relativity there have been many attempts to derive quasi-local conservation 
laws These conservation laws, if they exist, would serve as useful selection rules con- 

straining the future development of a spacetime, given the fact that the direct integration 
of the Einstein's evolution equations is practically impossible. Recall that in the Newtonian 
theory, the conservation of momentum J2p = const, immediately follows from Newton's 
third law, which is no more than the consistency condition implementing the second law. In 
general relativity, the consistency conditions for the evolution are already incorporated into 
the Einstein's equations via the constraint equations. Therefore it has been long thought 
that conservation laws should follow directly from the Einstein's equations [^,|^. The purpose 
of this paper is to show that this is indeed the case, and using the canonical (2,2)-formalism 
of the vacuum general relativity, I shall derive from the Einstein's equations on the time- 
like boundary surface the quasi-local energy conservation law, which relates the quasi-local 
gravitational energy loss to the net gravitational energy flux across a finite region of a given 
spacetime. 

Let us start from the following line element 

ds'^ = -2dudv - 2hdu^ + e^^p^b (^Z/" + A^du + A^dv) {dy^ + A^du + A^dv) , (1) 

where +, — stands for -u, w , respectively, and pab{ci, 6 = 2,3) is the conformal 2-metric of the 
transverse two-surface N2 defined hj u,v = constant, which satisfies the condition 

det pab = l. (2) 

Using the spacetime diffeomorphism invariance, one can always restrict the metric to the 
above form. Notice that these fields are functions of all the coordinates (m, since we 

assume no isometries. The integral Jq of the scalar curvature of the metric (|lD is, 

Iq = yiiu (if ?/ Lq + surface integral, (3) 

where Lq is given by 

Lo = —e^^'pabF^^F^' + e''{D+a){D^a) - ^eVV'(^+Pac)(l^-P6d) 

+e'^i?2 - 2e''{D_h){D_a) - he''{D_af + i/ieVV'(^-Pac)(/^-PM). (4) 
Here R2 is the scalar curvature of A'"2, and the notations are summarized below: 

F^^ = d+A_^ - d-A_^ - [A+, A_]l 

= d+A^ - d^A^ - A^dcA^ + A_'dcA_^, (5) 
D±a = d±(r - [A±,a]L 

= d±a - A^daa - daA£, (6) 
D±h = d±h-[A±,h]i^ 

= d±h - A^dah, 
D±Pab = d±Pab - [A±,p]-Lab 

= d±Pab - A^dcpab - {daADpcb - (ObA^) Pac + (<9c^±)Pafe, (7) 
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where [A±, *] is the Lie derivative of * along the vector field A± := A^da- The integral /q 
appears as an action integral of a Yang-Mills type theory defined on a fibre bundle, which 
consists of a (l+l)-dimensional base manifold coordinatized by {u, v), and the 2- dimensional 
vertical space N2 coordinatized by y"". In this picture, the Yang-Mills fields are valued in 
the Lie algebra associated with the group of the diffeomorphisms of N2 

In addition to the eight equations of motion that follow from (^) by variations, there are 
two equations associated with gauge- fixing the spacetime metric to the form (|^). These equa- 
tions, from which the quasi-local energy conservation law is shown to follow, are obtained 
by varying the Einstein-Hilbert action before the gauge fixing condition (|ip is introduced 
IPJ^. They are found to be 

(i) e"D+D^o + e^'D^D+a + 2e'"(D+a)(D_a) - 2e" {D ^h){D ^o) 
~e'y,,F^^F^^ - e-R2 - he-{{D_ay - ^pV'(^-Pac)(/^-P6d)} = 0, (8) 

(ii) -e^'Dla - -e^'iD+af - e''{D_h){D+a) + {D+h){D_a) 
+2he''{D_h){D_a) + e'^F^^dah - ^e>V'(^+Pac)(I^+P6d) + da{p'''dbh) 
+/ie-{ - + ipV'(^+Pac)(/^-Pfed) + ^e"^ PabF^^F^' + R^} 
+/i2e'^{(D_a)2 - lpV'(^-Pac)P-P6d)} = 0. (9) 

The above equations, being first-order in should be regarded as constraint equations 
associated with the partial gauge-fixing of the spacetime diffeomorphisms. Thus, in this 
formalism, the natural vector field with respect to which the evolution is to be defined is 
D^. Let us define the momenta ttj = {nh, Ha, T^a, tt"'^} conjugate to = {h, a, A_^, pab} as 

-I ■■= -Qrr, (10) 

where := D_q^ . Then we have 

Tih = -2e'^(Z}_a), 

TT^ = -2e'^(D_/i) - 2he"{D_a) + e'"(L)+a), 
vTa = e^^PabF+t 

vr- = he"p'''p'\D_pM) - ^e>V'(^+P6d), (H) 
where ■n"'^ is traceless since detpafe = 1. The Hamiltonian density Hq defined as 

■= TTjq' - Lo (12) 

is given by 

+l.n'^\D^Pa,) + i-e>'^V''(^+Pac)(/^+P6d) - e-R2. (13) 



3 



In terms of these variables the two equations (§) and become 



(i) Ho - d+TiH + da[A^'Kh + e-V'vTfc) = 0, (14) 
ii) -KhD+h + n^D+a + ix'^^'D+pab - d+ {2h vr;, + 2e''D+a) + da {2h ix^A^ 

+2A/e"D+a + 2he-''p''^Tn, + 2p"^9fc/i) = 0. (15) 



It is from the equation (15) that the quasi- local energy conservation law follows. To see 
this, let us consider the asymptotically flat spacetimes which approach the Schwarzschild 
solution 

ds^ — > -2dudv - f 1 - — ) du^ + v'^dVL^ (16) 

in the asymptotic limit. Here d/du is the timelike Killing vector field. Let us suppose that 
the spacetime we are considering is asymptotically fiat, so that N2 is a compact two-surface 
without boundary which is asymptotic to 5*2, and that the u-coordinate is asymptotic to the 
Killing time (see Figure 1). From the metric (|l|), we find that the hypersurface defined by 
(iw = is a null hypersurface, and the hypersurface defined by (if = is either a timelike or a 
spacelike hypersurface. Let us restrict our discussion to the region where = is a timelike 
hypersurface. The integral of the equation (|T5|) over a timelike hypersurface N2 x [mq, mi] 
becomes, with a suitable normalization. 



/ d^y (h -Kh + e^D+o) - i d^y (h txh + e^D+a) 

= [luh^y UhD+h + Ti^D+a + tt'^^D+p^^). (17) 

In general, the energy flux across a surface is given by the energy momentum tensor Toj, 
which is of the form 

Toi ~ T^^di<^ (18) 

in scalar field theories, where vr^ is the conjugate momentum of the scalar field 0. The 
integrand of the r.h.s. of (|1^) is precisely of this form, and therefore it represents the 
gravitational energy flux across the two-surface A^2, integrated from mq to u\. Then the 
l.h.s. of ( p!7D should represent the net change of the gravitational energy of the region 
enclosed by the compact two-surface N^. Thus our proposal of the quasi- local gravitational 
energy E{u^ v) of a region enclosed by a two-surface A^2 on the timelike surface defined by 
V = constant at the instant u is 

E{u,v) := i d^y(h7rh + e''D+a)+Eo{v), (19) 
iovr Ju ^ ' 

where Eq[v) is an undetermined substraction term which is also a two-surface integral. The 
equation ([T7|) then becomes 

E[u^,v) - E[uo,v) = flui d^y (nhD+h + 7T,D+(T + n'^'D+pab). (20) 

d/TT Juo Ju ^ ' 
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The flux integral in general does not have a deflnite sign, since it includes the flux carried 
by the in-coming as well as the out-going gravitational waves. But in the asymptotically flat 
region, the flux integral turns out to be negative-deflnite, representing the physical situation 
that there is no in-coming flux coming from the inflnity. 

For the asymptotically flat spacetimes, let us show that the relation (EOf) reduces to the 



standard Bondi mass-flux relation |jTO[ in the asymptotic limit. If we choose the substraction 
term Eq{oo) at inflnity as 

Eo(oo) = lim — hn (v^d^a), (21) 
the total energy E{u, oo) becomes 

E{u, oo) = m{u), (22) 

where m{u) is the Bondi mass of the asymptotically flat radiating spacetime. In order to 
show that the r.h.s. of (^) becomes the corresponding Bondi flux integral, we need to know 
the asymptotic fall-off rates of the flelds in a more detail. In the asymptotic region where 
an out-going twist-free null vector fleld exists, the metric is approximated as 

ds'^ — ^ -2dudv - 2hdu^ + e^pab [dy" + A^du) (dy^ + A_^du) , (23) 

where d/dv is the twist-free null vector fleld. If we compare this line element with the 
Schwarzschild spacetime (|16|), then the fall-off rates of the flelds are found to be 



f ^1 smt 



^){l + ^^^^^ + 0(l/.^)|, 

I 1J i 



p^^ = (sin 



Pi)^ = h 0{1 V ), 

V 

2/1=1-^ + 0(1^), 



A^^ = Oil/v), (24) 

where a, (3, 7, and 5 are arbitrary functions of {u, ^, if). Using these conditions, the r.h.s. of 
(pop becomes. 



62tc Juo Ju ^ 

= flucfdM^ I - ^eVV'('9+Pac)(9+pbd) + e'^(9+a)^ - 2he''{d+a){d^a)}. (25) 

OZTT Juo Ju ^ 2 J 

As is shown below, the integral (|25| ) can be made negative-deflnite, if we use the residual 
symmetry of the metric (p3|) which consists of, apart from the obvious diffeomorphisms 
of N2, the following coordinate transformation 
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V — > v' = V + e{u,^,(p), (26) 

while keeping {u,-&,(f) constant. Here e is an ar^rfrary function of {u,d,(p). This symmetry 
transformation reflects the freedom of choosing the origin of the afline parameter v in an 
arbitrary way. In the new coordinate system, we find from (|24D that the asymptotic form of 
e°" becomes 



^(sin^9){l + i(a - 2e) + 0{l/v'^)]. (27) 
Let us choose the function e 



such that becomes 



Then we have 



6 = ^a, (28) 



;sin?9){l + 0(l/i;'2)}. (29) 



d+a = 0{l/v'^), d-a = 0(l/f'), (30) 
so that the last two terms in the r.h.s of ( pSj ) vanish in the limit v' —* oo. Therefore the 



energy difference (|20|) becomes 

1 /""i f 

E{ui,oo) - E{uo,oo) = -— / duidnv^p'''>p'\d+pac){d+pM) < 0. (31) 

D47r Juo Ju 

This relation is precisely the Bond mass-loss formula, relating the mass loss in terms of a 
negative-definite flux of the out-going gravitational waves. The r.h.s. of the integral is a 
bilinear of the gravitational current j\ 

j\ = p'^'d+p,,, (32) 

which is precisely the shear degrees of freedom of the gravitational field. 

Let us now consider a spherical ball of radius v, filled with a perfect fluid with energy 
density p{v), and find the total energy contained in the ball. The Einstein's equations are 
modified by the presence of the fluid, and the solution is given by 

ds^ = -2dudv - 2h{y)du^ + v'^d^^, (33) 

where 

2h{v) = 1 - m{v) = ATcTdv'v'^piv'). (34) 



The substraction term at the infinity Eo{oo) in (^) suggests that the substraction term 
Eq{v) at a finite distance should be chosen as 

Eo{v) = -^<fdn(v''d^a). (35) 

lOTT JS2 ^ ' 
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Then E{u, v) becomes 

E(u,v) = — (fdQ(l-2h)( v^d^a 
167r Jso ^ ^ ^ 



S2 

= m(f), (36) 

which is identical to the Misner-Sharp energy for spherically symmetric distribution of 
an ideal fluid. However, this is not the proper mass. The proper volume element dV of the 
u = constant null hyper surf ace is given by 

{2h)-^/\^dvdQ, (37) 

so that the proper quasi-local mass mp{v) is given by 







4Wrf.V^(l-^)-%(.'). (38) 



/o ^ v 

The difference between E{u,v) and mp{v) can be interpreted as the gravitational binding 
energy Eb{v), so that E{u,v) becomes 

Eiu,v) = mp{v)-EB{v). (39) 

In the Newtonian limit where E{u,v)/v is small, Eb{v) becomes 

E^{v) ^ 4nrdv'v''p{v')(^^^), (40) 

which is just the total Newtonian binding energy of mp{v') interacting with the fluid shell 
of mass indv'v''^ p{v'), integrated over the ball. Notice that Eb{v) has the correct sign [|I|, 
so E{u,v) indeed represents the total gravitational energy within the spherical ball of the 
radius v, including the binding energy. 



It is worth mentioning that from the equation (|T^) the quasi-local conservation of mo- 
mentum also follows. The resulting equation has a similar structure with the Euler equation 
in fluid dynamics []T2[. The Hamiltonian density Hq in this (2,2)-formalism, which corre- 



sponds to the momentum flux, consists of terms which are at most quadratic in ttj. From this 
point of view, terms linear in tt/ may be regarded as viscosity terms, and terms independent 
of TT/ as pressure terms. 

In summary I derived a quasi-local version of the mass-flux relation for the vacuum gen- 
eral relativity from the Einstein's equation on the timelike boundary hypersurface. The 
proposed quasi-local mass and flux integral are shown to reproduce the standard results in 
the asymptotic region, and agrees with the Misner-Sharp energy for the spherically sym- 
metric fluid in a finite region. Unlike other proposals of quasi-local gravitational energies, 
our expression is canonical, being derived from the first-order (2,2)-formalism. From the 
viewpoint of the correspondence principle, our derivation is particularly satisfactory, since 
the quasi-local conservation in this formalism comes from the "consistency equation" , as in 
the Newtonian theory. 
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95-0702-04-01-3. 
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FIGURES 




FIG. 1. A region of spacetimc where the surfaee = is a null hypersurface, and the surface 
= is a timelike hypersurface. Their intersection defines A^2- The vector field d/du is normal 
to if and only if A,^ = 0, and its norm approaches —1 near the infinity. 



9 



